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Abstract 

In this paper we analyze the structure of supersymmetric vacua in compactifications of the 
heterotic string on certain manifolds with SU(3) structure. We first study the effective théories 
obtaincd from compactifications on half-flat manifolds and show that solutions which stabilise 
the moduli at acceptable values are hard to find. We thcn dérive the effective théories associated 
with compactification on generalised half-flat manifolds. It is shown that thèse effective théories 
are consistent with four-dimcnsional N = 1 supergravity and that the supcrpotcntial can be 
obtained by a Gukov-Vafa-Wittcn type formula. Within thèse generalised models, we find 
consistent supersymmetric (AdS) vacua at weak gauge coupling, provided we allow for gênerai 
internai gauge bundles. In simple cases we perform a counting of such vacua and find that a 
fraction of about 1/1000 leads to a gauge coupling consistent with gauge unification. 
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1 Introduction 



There is now a considérable body of work on moduli stabilization, facilitated by flux of the Neveu 
Schwarz-Neveu Schwarz (NSNS) and Ramond-Ramond (RR) anti-symmetric tensor fields, in the 
context of type II théories. Specifically, within type IIB it has been shown that a combination of 
NSNS and RR three-form flux can stabilize the dilaton and ail complex structure moduli, while the 
Kàhler moduli have to be fixed by other effects such as non-perturbative contributions [2] or perhaps 
higher-order a' corrections [3J. The consistency of thèse procédures, including the interplay between 
a' and non-perturbative corrections, was analysed in Refs HJ |H1 El El Qi3 Within type IIA 
théories, on the other hand, both odd and even degree form field strengths are available, so that flux 
potentials for complex structure moduli as well as Kàhler moduli will typically be generated [111 112] 
(see also Réf. |13j for N = 1 models). One may therefore hope that ail moduli can be stabilised 
by flux in some such models and spécifie examples have indeed been found |14l 1151 1161 1171 118j . 
although it appears that in generic models of this kind some flat directions are still left over. 

Traditionally, the heterotic string has been considered the most attractive string theory, with 
the présence of (preferably Eg x E$) ten-dimensional gauge fields leading to a large number of 
supersymmetric compactifications with phenomenologically interesting properties [HJ- It bas also 
been known for a long time that heterotic NSNS three-form flux can stabilize ail complex structure 
moduli of the theory [201 [^. More recently, this subject was addressed in Refs [2211221. However, 
in the absence of any further (RR) antisymmetric tensor fields, the potential for stabilizing the 
remaining moduli seems rather limited compared to type II théories. This apparent problem can 
be overcome by departing from Calabi-Yau compactifications and by considering the heterotic 
string on gênerai manifolds with SU(3) structure. Such models were analysed in Refs j24l 1251 126[ 
1271 1281 1291 I3U| where gênerai aspects of compactifications on non-Kàhler manifolds were studied. 
Recently, a more gênerai analysis which takes into account the effects of a gaugino condensate 
in ten dimensions appeared in Réf. [35- The generic form of the superpotential was inferred in 
[231 Oïïl l2"7j , but its detailed analysis was not possible due to the lack of knowledge of the moduli 
space of thèse manifolds. On the other hand the low energy effective action and an explicit form 
for the superpotential in terms of the low energy fields were found in Réf. [32], where half-flat 
mirror manifolds were used as compactification spaces. Such manifolds arise in the context of 
type II mirror symmetry with NS fluxes |33| and, in some appropriate région, their moduli space 
was conjectured to be similar to that of a normal Calabi-Yau manifold. This conjecture has been 
applied in Réf. [221 to dérive the low-energy theory for thèse compactifications and, in particular, 
the superpotential as an explicit function of the moduli fields. In particular, it was found that 
the intrinsic torsion of the half-flat mirror manifolds gives rise to a superpotential for the Kàhler 
moduli. Thèse results suggest that, by combining the intrinsic torsion of sufficiently gênerai classes 
of SU(3) structure manifolds with NSNS flux, moduli stabilization in heterotic compactifications 
can be as flexible as in type II models. 

This is precisely the line of work we would like to further develop in the présent paper. We 
will first study in détail supersymmetric moduli stabilization in the heterotic string on half-flat 
mirror manifolds, based on the effective théories of Réf. [H21- As we will see, the torsion of half-flat 
manifolds and the allowed H-Ruxes are insufficient to fix ail Kàhler and complex structure moduli. 
We therefore move on to the more gênerai class of manifolds with SU(3) structure described in 
Refs |341 1351 131)] . which we will refer to as generalised half-flat manifolds. For this class of spaces we 
first show, by an explicit réduction of the bosonic action, that the heterotic Gukov-Vafa-Witten type 
formula for heterotic compactifications [321 1371135] leads to the correct resuit for the superpotential. 
This resuit is then applied to a detailed analysis of moduli stabilization for those models. 

The flux and torsion superpotential, W, for ail models considered in this paper is a function 
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of the Kàhler moduli T % and complex structure moduli Z a , but it turns out to be independent of 
the dilaton S. Hence, the dilaton is not stabilised at this stage. However, in the Es x Es case, one 
expects hidden sector gaugino condensation to generate a non-perturbative dilaton superpotential 
which should be added to W. We will, therefore, use this non-perturbative contribution to stabilize 
the dilaton. It turns out that, in order to fix S at sufficiently weak coupling (and to be in the 
large radius and large complex structure limits) we need global minima for the Kàhler and complex 
structure moduli which correspond to superpotential values Wo with |Wq| <C 1. This is quite 
analogous to a similar requirement in type IIB models ]2], where it is necessary to ensure moduli 
stabilization at large radius. The original models of heterotic gaugino condensation with flux [201121] 
were discarded precisely because this condition was difficult to satisfy due to the quantization of 
fluxes. However, we find that cancellations leading to small |Wo| are possible for our generalised 
models. We carry out a statistical analysis in those cases, counting the number of vacua as a 
function of |Wo| and the maximal flux value. As the value of |Wo| détermines the value of the 
dilaton, this counting analysis is directly relevant to the question of how many vacua realize a 
phenomenologically acceptable gauge coupling. 

The outline of the paper is as follows. In Section |21 we briefly review the low energy effective 
theory of the heterotic string on half-flat mirror manifolds [35]. In addition, we work out the 
generalization of this effective theory expected for the more gênerai spaces proposed in |34j . We 
will show explicitly that the potential obtained from compactification (which includes a part from 
the non-vanishing scalar curvature of the internai space) can be obtained from a Gukov-Vafa-Witten 
type superpotential for manifolds with SU(3) structure, which was derived in [32]. In Section UJ 
we set up our four-dimensional models in a way suitable for the discussion of moduli stabilization 
which includes gaugino condensation and flux quantization. This section is largely self-contained 
and the reader mostly interested in the four-dimensional aspects of our analysis may want to 
skip Section [21 and move on to Section |3J straight away. Moduli stabilization within models based 
on half-flat mirror manifolds [32] is discussed in Section 0] In Section [5] we discuss the models 
based on the more gênerai half-flat spaces introduced in Section |2J We conclude in Section El 
Various technical détails are deferred to the three appendices. Appendix[Â]contains a calculation 
of the scalar curvature of the generalised half-flat spaces, which is essential in establishing the 
consistency of the generalised models of Section [2] In Appendix[B] we have collected a number of 
useful relations on spécial geometry, while Appendix O summarizes our four-dimensional N = 1 
supergravity conventions. It also includes an elementary proof that supersymmetric AdS vacua of 
this theory are always stable. 

2 The heterotic string on half-flat manifolds 

In this section we will review the compactification of the heterotic string on half-flat mirror mani- 
folds jH2] and présent an extension of this work to the spaces proposed in Réf. [BU EE1 Œ] • 

2.1 The heterotic string on half-flat mirror manifolds 

Half-flat mirror manifolds arise in the context of type II mirror symmetry |33j and can be thought 
of as mirror duals to Calabi-Yau manifolds with NSNS flux. More specifically, given a mirror pair 
(X, Y) of Calabi-Yau manifolds, the mirror of, say, IIB on Y in the présence of NSNS flux is IIA on 
a half-flat mirror manifold X (without flux). This half-flat mirror manifold X is closely related to 
the original Calabi-Yau mirror X in that it can be characterized by the two Hodge numbers h^ 1 ' 1 ' 
and M 2 ' 1 ) of X and carries sets of two- three- and four-forms analogous to the sets of harmonie 
forms on the Calabi-Yau space X. Specifically, on X we dénote by (uii) and {Cb 3 ) a basis for the 
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two- and four-forms respectively, where ■ ■ ■ = 1, . . . , h^ l,1 \ which satisfy 

u)iMV=6{. (2.1) 



x 

Further, on X, one can define a set of symplectic three forms (et a, P B ) where A,B, - ■ ■ = 0, . . . , hS 2,1 ' 
with 

/ a A A{3 B = ÔA, [ et a Aqb = / (3 A A (3 B = . (2.2) 

Being manifolds with SU(3) structure [25], half-flat mirror manifolds carry a two-form J and three- 
form Q which are the analog of the Kàhler forai and the holomorphic (3, 0) form on Calabi-Yau 
manifolds . As on Calabi-Yau manifolds thèse forms can be expanded as 

J = Put , (2.3) 
fl = Z A a a -Ç A f3 A , (2.4) 

where t % and Z A are the équivalent of Kàhler and complex structure moduli. As usual, the coefficient 
Qa can be obtained from a holomorphic pre-potential Q = Q(Z ), homogeneous of degree two, as 

Sx = . (2.5) 

So far, the set-up has been exactly as for Calabi-Yau manifolds. The main différence is that the 
forms (u>i) and (otA,P A ) are no longer harmonie but rather satisfy 

doJi = eiP° , da = , da a = df3 A = , dû i = . (2.6) 

Here are h^ 1 ' 1 ^ parameters (real numbers) which characterize the torsion of the half-flat mirror 
manifold under considération. 

Having described the basic structure of half-flat mirror manifolds, let us now review the com- 
pactification of the heterotic string (at lowest order in a') on those spaces. Besides the metric, 
there are two other bosonic fields, namely the dilaton s = exp(— 2(j>) and the NSNS two-form B. 
The latter can be expanded as 

Ê = B + r'ui , (2.7) 

where B is a four-dimensional two-form which can be dualised to a scalar a and r* is a set of axions. 
Together with the dilaton and the Kàhler moduli, thèse fields pair up into four-dimensional chiral 
multiplets as 

T = t 1 + it , (2.8) 
S = a + is. (2.9) 

In ternis of the projective coordinates Z A the complex structure chiral multiplets Z a , where 
a, b, ■ ■ ■ = 1, . . . , M 2,1 ), are obtained by Z a = Z a /Z° and we write thèse fields as 

Z a = C + iz a - (2.10) 

The Kàhler potential for those fields in the large radius limit is then of the standard Calabi-Yau 
form, that is, 

K = K^+K^ + K^ , (2.11) 

1 Although the manifolds discussed in this paper are generally neither complex nor Kàhler, we will frequently use 
Calabi-Yau terminology and, for example, refer to J as Kàhler form. 
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with 



= -ln('|jC 



K (S) = -i n (i(S-S)) , 



and 



K = -d iih {T i -f i ){T> -P)(T k - f k ) = d ijk tHH k , (2.13) 

o 

£ = j{z A g A -z A g A ) . (2.14) 

Here, d^-fc are numbers analogous to the intersection numbers of the associated Calabi-Yau space 
X. Later, we will be working in the large complex structure limit, where the pre-potential G can 
be written as 

1 d ahc Z a Z b Z c 

Q = — 6 Yo , (2-15) 

with d a bc analogous to the intersection of the associated mirror Calabi-Yau space Y. In this case, 
the complex structure Kàhler potential takes a form similar to the one for the Kàhler moduli, that 

K = l -d abc (Z a - Z a )(Z b - Z b ){Z c - Z c ) = d abc z a z b z c . (2.16) 



Let us now discuss the superpotential. In Réf. |32| it has been shown that, for gênerai heterotic 
compactifications on manifolds with SU(3) structure, the superpotential to order a' can be obtained 
from the Gukov-Vafa-Witten type formula 

W = V8 QA{H + idJ), (2.17) 
Jx 

where H is the NSNS field strength. For half-flat mirror manifolds this field strength can be written 

as 

H = dB + dr* A oJi + T^iP + F flux . (2.18) 

where the first three terms have been computed by taking the exterior derivative of Eq. (|2.7|) . Note 
that the third term is new and arises because the forms uji are no longer closed, see Eq. (|2.6|) . We 
have also added on an additional NSNS flux contribution 



#flux = (/i a a« - e a f3 a ) , (2.19) 

with electric and magnetic flux parameters e a and fi a , respectively. If we arrange the a' terms 
in the H Bianchi identity to cancel (for example by choosing the standard embedding) if follows 
that dH = 0. For this reason we have dropped the term proportional to the non-closed form aç, 
in Eq. I|2.19|) . We have also omitted a possible term proportional to in (|2.19f) which can be 
absorbed into a re-definition of the axions r', as is évident from Eq. (|2.18|) . Inserting the field 
strength 1)2.18(1 . the (3, 0) form 1)2.4)1 and dJ = e^P into the gênerai formula 1)2.17)1 . one finds the 
superpotential 

W = VsietT + e a Z a - ^ a Ç a ) , (2.20) 
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where the basic intégrais (|2.2|) have been used. This resuit has been checked in Réf. where the 
four-dimensional scalar potential was calculated from an explicit réduction of the ten-dimensional 
bosonic action of the heterotic string. This scalar potential has three contributions which arise from 
the third term in Eq. (|2.18|) and the NSNS flux (|2.19|) . both inserted into the H kinetic term, and 
the non-vanishing scalar curvature of the half-flat mirror manifolds. Thèse three contributions lead 
to a potential which can be exactly reproduced from the above superpotential, using the standard 
relations of four-dimensional N = 1 supergravity (see Appendix [U] for a summary of supergravity 
conventions). In the following subsection we will generalize this calculation to a larger class of 
manifolds with SU(3) structure. 

2.2 Setup for the extended models 

Having discussed the basic models obtained from the compactification on half-flat mirror manifolds 
we can now study a généralisation of the half-flat spaces which was proposed in Réf. [21]. The same 
class of spaces appeared in [23 EU and it was argued to be the correct Ansatz for a consistent 
Kaluza-Klein truncation to four dimensions. Here we will use the prescriptions given in the above 
références, and show that such a truncation is indeed consistent with supersymmetry. In particular, 
we will show that the four-dimensional scalar field potential derived from a compactification on 
those generalised spaces is consistent with the Gukov-Vafa-Witten type formula (|2.17|) for the 
superpotential. 

We start by reviewing the main features of this new class of manifolds with SU(3) structure 
which we will dénote generalised half-flat manifolds. We will mostly follow Refs. |34[ I36j . As we 
have done for the half-flat mirror manifolds, the existence of two-forms (wj), four-forms and 
three-forms (ctA,P B ) satisfying the basic intégral relations 1)2.1(1 and (|2.2() is postulated. However, 
the crucial differential relations (|2.6fl are now generalised to 

dtui = p Ai (3 A - qfa A , da A = VAi& , d(3 A = qftf , dw* = , (2.21) 

with (real) torsion parameters pAi and qf. From d 2 Wi = one concludes that the additional 
constraints 

PAiqf ~ qfpAj = , (2.22) 

have to be imposed, for a consistent définition of the exterior derivative. In what follows the 
above defining relations are enough in order to dérive the low energy action which arises from 
compactification on thèse manifolds. In appendix^we will have more to say about the geometry 
of thèse spaces and, in particular, about their torsion classes, which differ from those of a half-flat 
manifold. 

The expansion of the Kàhler form, J, the (3,0) form, £1, and the NSNS two-form, B, in terms 
of the basic forms remains unchanged and is given in Eqs ((2.3|) . (|2.4j) and (|2.7|) . This also means 
that we have the same set of moduli fields, 2 namely the Kàhler and complex structure moduli T % 
and Z a and the dilaton S. Whenever exterior derivatives are taken we now have to work with the 
generalised relations (|2.21|) . This means that the NSNS three- form field strength associated to (|2.7|) 
is given by 

H = dB + dr i A + Ti(p Ai p A - qfa a ) + # flux , (2.23) 

2 Strictly speaking the fields we are talking about are no longer moduli as the potentials generated are not flat in 
thèse directions. However we continue to call them moduli in order to stress that we are interested in the fields which 
were the moduli of the related Calabi-Yau compactifications. 
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where, as before, we have added on the NSNS flux part 



flux 



H a A ~ e A /3 J 



(2.24) 



with electric and magnetic flux parameters e A and \x A . If the RHS of the heterotic Bianchi identity 



o 



dH= — (tr(F A F) - tr(R A R)) , 



(2.25) 



vanishes (for example, by choosing the standard embedding), then H needs to be closed which 
implies the further constraints 

Ai. - CAqf = , (2.26) 

between flux and torsion parameters. On the other hand, the RHS of Eq. (|2.25|) . although neces- 
sarily exact, can be non-zero, so that the constraint ()2.26|) can be avoided by, for example, more 
complicated choices of the gauge bundle. It is convenient to introduce the following combinations 



îa = eA - T l p Ai 



(2.27) 



of fluxes, torsion parameters and Kàhler moduli in terms of which the NSNS field strength can be 
expressed as 

H = dB + dT l u t + Re(fi A )a A - Re(ë A )(3 A . (2.28) 
For the exterior derivative of the Kàhler form J one finds 



dJ = tfduji = lm(fl A )a A - lm(ë A )(3 A , 



(2.29) 



where the differential relations (|2.21|) and the définitions (|2.27|) have been used. Thèse last two 
results for H and dJ, together with the standard expansion for the (3, 0) form (|2.4j) and the 



basic intégrais (|2.2j) . can be used to evaluate the formula (|2.17j) for the superpotential. A simple 
calculation leads to 

W = V8(ë A Z A - £l a Ç a ) . (2.30) 
We will now verify this resuit by an explicit réduction of the ten-dimensional bosonic action. 



2.3 Réduction for the generalised models 

The starting point for the compactification is the lowest order in a' of the bosonic part of the 
ten-dimensional effective action of the heterotic string. This is given by 



Si 



1 



0, bosonic 



2^10 Jm w 



-2d> 



R*l- 4d(f> A *dc[> + -H A *H 



(2.31) 



As the main assumption for compactifications on generalised half-flat manifolds is that the light 
spectrum of normal Calabi-Yau (and also half-flat) compactifications is unchanged, we will not 
be concerned with the dérivation of the kinetic terms for the various fields one obtains in four 
dimensions. They are exactly as discussed for the case of half-flat mirror manifolds, see Eqs (|2.11[) - 
(|2.14|) . Instead we concentrate on the scalar potential. As explained in the previous section, one 
contribution to the four-dimensional potential arises from the H kinetic term with (|2.28j) inserted. 
A standard calculation jlll I4(J| leads to 



e~ K V H = -4e~ K(Z) [Re(ë A ) - Re{f)M C A] (imAT 1 )^ [Re(ê A ) - Re(fi c )M C A] • (2.32) 



\\AB | 
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Here the matrix Ai is the period matrix ()B.7|) which, for the complex structure sector, is also given 
by the relations (jB. 10|) . 



The second contribution arises from the Einstein Hilbert term in ()2.31|) and is due to the non- 
vanishing scalar curvature of the half-flat spaces. The calculation of this scalar curvature, for the 
spaces characterized by the relations (|2.21|) . is somewhat non-trivial and has been carried out in 
AppendixEI The resuit is 

e- K V R = -4e- K(Z) [lm(6 A )-Im(/i c )Al CA ](lmAl- 1 ) AB [lm(è A )-Im( / i c )^^ 

+8E i Ë j {g ij -4tH j ) , (2.33) 

where we have introduced the notation 

Et = PAi Z A - q A Q A . (2.34) 

It is not hard to see that, provided the constraints (|2.22|) and (|2.2(i|) are satisfied, the total potential 
takes the form 

e~ K V = -4e- K(Z) [è A - ï c M C a] (imAT 1 ) AB [è A - fi c M CA ] + S^Ë^ - . (2.35) 

We now need to verify that this potential indeed originates from the superpotential (|2.3U|) via the 
standard supergravity formula (|C.2j) . Since the index X in this formula runs over ail chiral fields 
which, in our case, consist of the dilaton S, the complex structure moduli Z a and the Kàhler moduli 
T % , we will discuss each case separately. 

First of ail notice that, since the superpotential (|2.3U|) does not dépend on S, the contribution 
of the dilaton-axion chiral superfield to the potential can be found from (|2.12|) to be simply 

e~ K V s = \W\ 2 . (2.36) 



For the complex structure moduli we obtain 

-j=D a W = (è B - Ç B ciï C )D a Z B , (2.37) 

where we define Gbc = dsdcQ- Using the relations ()B.9|) and ()B.6|) one immediately finds 

g al D a WD- b W = -Ae~ K(Z) (imAi" 1 )^ (è A - jf ' M CA ) (ë B ~ fi D M DB ) - 8\è A Z A - fi A Q A \ 2 . 

(2.38) 

Note that the first term in the above équation is similar to the first term of (|2.35|) . except for the 
complex conjugations which do not work out quite right. However, it is just a matter of algebra 
to show that thèse two terms are indeed identical provided that the constraints (|2.22|) and ()2.26j) 
hold. 

Also, the second term in (|2.38|) very much resembles the square of the superpotential, but here 
the complex conjugations can not be exchanged so easily. In turn one obtains 

g al D a WD- b W = -4e- /<(Z) (lm^- 1 ) AB (è A - f2 c M CA ) {f B - Ï D M DB ) - \W\ 2 + Y , (2.39) 
where by Y we have denoted the combination 

Y = 32/ '/'/•;,/•.) - 2V8it i (E i W - E { W) . (2.40) 



7 



Let us finally deal with the Kàhler moduli contribution to the N = 1 potential. Using formulae 
(|B.15|) - (|B.19|) on the Kàhler moduli space we find 



g ij DiWDjW = g ij diWdjW + 3|iy| 2 + 2U i (Wd î W - WdiW) . (2.41) 

To this end it is useful to make the dependence of the superpotential (|2.3U|) on the Kàhler moduli 
explicit by writing 

W = V%{-E t T l + e A Z A - n A g A ) , (2.42) 
where Ei were defined in 1)2.34(1 . Hence we have 

d t W = -V8E { . (2.43) 

With this we see that the last terms in Eqs (|2.41() and (12.40(1 cancel identically. Moreover, the \ W\ 2 
terms from équations 1)2. 36 j) . (|2.39j) and 1)2.41)1 cancel against the — 3|VF| 2 in équation (|C2|) while 
the remaining terms precisely combine into (|2.35() . This concludes our dérivation of the potential 
1)2.35)1 from the superpotential 1)2. 30 J) , and establishes a strong argument for the consistency of the 
compactifications on the manifolds presented in section l2~2*l which were introduced in Refs. |34l 1351 

nu. 

We conclude this section by comparing the superpotential (|2.3U|) which we have just derived 
with the one obtained in type IIB compactifications. There, the fluxes are "complexified" in a way 
that involves the IIB complex coupling. In our case, the flux parameters are "complexified" to è 
and fl in Eq. (|2.27|) due to their dependence on the Kàhler moduli. Apart from this "exchange" 
of Kàhler moduli and dilaton, the resemblance between the two superpotentials is quite striking. 
This confirms our expectation that heterotic théories can be as flexible with regard to moduli 
stabilization as type II théories when non-trivial torsion is included. 



3 General structure of low-energy théories 

So far we have concentrated on how four-dimensional models arise from compactifications of the 
underlying ten-dimensional theory. In the remainder of the paper we will analyze the implications 
of thèse four-dimensional models for moduli stabilization, and the purpose of this section is to set 
up ail the necessary ingrédients, in a way that is convenient for this analysis. 

3.1 The models 

From now on, we will adopt the "phenomenological" définition of the chiral superfields in terms 
of its components, where the real parts are the "geometrical" moduli and the imaginary parts the 
axions. With respect to our previous convention, this corresponds to the simple transformation 
(j) X — ► —i(j) X (together with a sign flip of the axions) on ail fields. Explicitly, this means we are 
replacing the field définitions (|2~%|) . (f23|) and (|2~TU|) by 

S = s + ia , (3-1) 
r = f + iT l , (3.2) 
Z a = z a + i( a ■ (3.3) 

While our gênerai calculation for the four-dimensional effective theory was valid for ail values of 
the complex structure moduli we will, in the following, focus on the large complex structure limit. 



S 



This means that, from Eqs (|2.11|) — (|2.16|1 together with the above field re-definition, the Kàhler 
potential is given by 3 

K = — ln(S + §)- ln(8/C) - ln(8K.) , (3.4) 

with 

K = \d i]k {V + T){Ti + P)(T k + f k ) = d ijk t i tH k , (3.5) 

o 

IC = ]-d abc (Z a + Z a )(Z b + Z b )(Z c + Z c ) = d abc z a z b z c . (3.6) 

o 

Recall that dij k and d a b c correspond to the intersection numbers of the associated Calabi-Yau space 
and its mirror, respectively. Both the Kàhler and complex structure parts of the Kàhler potential 
are given in ternis of spécial geometry pre-potentials which, due to large radius and complex 
structure, are determined by cubic polynomials. The cubic nature of the pre-potentials means both 
moduli spaces constitute examples of very spécial geometry. Some useful relations for very spécial 
geometry, which we will apply subsequently, are collected in AppendixlBl 

Let us now turn to the superpotential. Inserting the derivatives (|2.5|) of the large-complex 
structure pre-potential (|2.15|) into Eq (|2.3()|) . along with the définitions (|2.27j) of the complex flux 
parameters, the explicit form of the superpotential W turns out to be 

W = -i(e -iTp 0i ) + (e a -iTp ai )Z a + l -(f, a -iT i qt)d abc Z b Z c 

+ -iTq° i )d abc Z a Z b Z c . (3.7) 
b 

As we have pointed out in section |2~21 the parameters in this superpotential are not independent 
but satisfy 

PAiqf ~ qfpAj = , (3.8) 
e A qf-v A PAi = 0. (3.9) 

Note that the first of thèse constraints follows from the property d 2 = of the exterior derivative 
and is, therefore, strictly necessary. The second one is a conséquence of dH = 0, which is the 
correct form of the heterotic Bianchi identity if the a' corrections on the RHS of Eq. (|2.25|) cancel 
by themselves, for example, by choosing the standard embedding. However this need not be the 
case, so that this second constraint can be avoided. 4 In this paper, we will study both cases with 
and without the second constraint. Finally, note that half-flat mirror manifolds correspond to the 
spécial case where we set eo = p a i = qf = M° = Qi = and poi = — &i m the superpotential ()3.7j) . 
This leads to 

W = ei T + e a Z a + % -d ahci x a Z b Z c , (3.10) 

which is the large complex structure limit of Eq. (|2.2U|) . as it should. 

The above Kàhler potential and superpotential feed into the gênerai formula for the four- 
dimensional N = 1 supergravity potential and we have summarized the relevant conventions in 

3 Several numerical factors from 12. lit to 12.161 were absorbed into the superpotential and the définition of the 
flux parameters in order to make the calculations in the following sections more straightforward. 

4 Since the calculation in section relies on équation 13.91 1 we may argue that this constraint cannot be relaxed. 
However, if we were to incorporate consistently ail the terms which appear at order a' , we would expect to find the 
same superpotential as before. In fact, this is precisely what the Gukov-Vafa-Witten formula 12.17H evaluated for a 
field strength H which includes the a' corrections predicts. 
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AppendixIO In this paper, we will only be concerned with supersymmetric vacua of thèse poten- 
tials, that is, solutions to the F-equations. Generically, such solutions have a négative cosmological 
constant (|C5|) and so they lead to four-dimensional AdS vacua. It is known [U that such vacua 
are always stable and Appendix[U]also contains an elementary proof of this fact. 



3.2 Gaugino condensation 

As the dilaton does not appear in the superpotential W, Eq. (|3.7|) . the potential will usually be 
runaway in this direction. Hence, if we want to have any chance of stabilizing ail moduli, we should 
consider additional contributions. As has been shown in Réf. [32]) the gauge kinetic function / of 
the four-dimensional gauge group SO(10) (g) E% for heterotic compactifications on half-flat mirror 
manifolds is given by 

f = S, (3.11) 

to leading order. Clearly this resuit extends to the generalised half-flat manifolds discussed in the 
previous section and more gênerai gauge bundles. Hence, hidden-sector gaugino condensation [20] 
leads to an additional dilaton-dependent superpotential term which is precisely what we need. We 
will, therefore consider the superpotential 

W = W + ke~ cS , (3.12) 

with W as given in Eq. (|3.7|) . Here k and c are constants, the latter being determined by the 
one-loop beta function of the gauge group. To make this more précise, we normalize the real part 
of the dilaton, s, such that 

1 4vr , 

s = = — , (3-13) 

a YM g Y M 

where çym is the Yang- Mills coupling constant. In terms of the one-loop beta function coefficient b, 
the constant c can then be written as c = 6ir/b. For gauge group E%, one finds b = 90 and, hence, 

= P.14) 

The pre-factor k is hard to fix precisely not least because corrections due to the two-loop beta- 
function will lead to an S'-dependent pre-factor of the exponent in W, which we neglect in the 
présent context. We will simply parameterize k as 

k = . (3-15) 
Va 

where a' is the string tension and k is a dimensionless constant which one expects to be of order 
one. 



3.3 Quantization of flux and torsion 

We would now like to be somewhat more spécifie about the quantization of flux and torsion param- 
eters in the superpotential. For the genuine fluxes this is easy to achieve |21j by imposing that H 
is an élément in the intégral cohomology (modulo normalization factors). It is less straightforward 
to see how the torsion parameters of the internai manifold should be quantized. For the half-flat 
mirror manifolds, this will be done via the mirror symmetry relation which was used in order to 
establish the existence of such spaces in the first place. Unfortunately, such a correspondence is not 
known for the more gênerai manifolds described in section 12.21 so we will have to make a plausible 
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assumption about quantization for thèse spaces, generalizing from the results obtained for half-flat 
mirror manifolds. 

Before we can find the quantization rules for the flux parameters we should fix the normalization 
of our moduli fields. We recall that the above models have been derived and are valid in the large 
radius and large complex structure limit. Hence, we adopt a normalization of fields where thèse 
limits correspond to field values 

tf > 1 , z a > 1 . (3.16) 

What does this convention imply for the underlying internai geometry? Recall that the dimension- 
less Kàhler moduli fields measure the volume of the various Calabi-Yau two-cycles C\ in units of 
some (six-dimensional) référence volume v. More precisely, we have 



.1/3 



J , (3.17) 



where J is the Calabi-Yau Kàhler form. In order to assure that t l > 1 indeed corresponds to the 
limit in which the "radius" of thèse cycles is bigger than one in string units, one has to fix this 
référence volume to be 5 

v = (4ir 2 a f . (3.18) 

In order to fix the normalization of the complex structure moduli in a similar way, it is useful 
to consider the mirror picture. The fields z a measure the size of two-cycles on the mirror and large 
radii for thèse two-cycles corresponds to large complex structure in the original model. The volume 
of thèse mirror two-cycles should be measured in units of the same référence volume ()3.18j) . in order 
for the two four-dimensional effective théories from the original model and its mirror to be identical 
(and the mirror map being trivial on the four-dimensional fields). With this convention, it is then 
clear that z a > 1 indeed corresponds to the large complex structure limit. 

We are now ready to discuss the quantization of the flux parameters which appear in the 
superpotential (|2.20[) , The basic quantization condition for the NSNS three-form field strength H 
is given by |2J 

/ H = 4vrV • l , (3.19) 
Jc 3 

where l is an integer and C3 represents any three-cycle in the integer homology. Following the 
explicit calculation of the four-dimensional potential by dimensional réduction in Réf. [32], it is 
easy to see that this quantization rule, together with the scale convention (|3.18j) . implies that 

e a = ^Lè a , (3.20) 
7rv a 

fx a = ^=ï\ (3.21) 



where e a and jl a are integers. Note that the counterintuitive numerical factors include the redéf- 
initions of the flux parameters, which were needed in order to rewrite the Kàhler potential and 
superpotential in the simpler form of (|3,4|) to (|3.7j) . In order to fix the quantization of the electric 
torsion parameters, e^, we should again consider mirror symmetry. On the mirror, thèse electric 
torsion parameters become electric flux parameters of the NSNS form. Given that our basic choice 



5 Of course, there is always an ambiguity of factors of 2ir in this calculation which cannot be easily fixed. To arrive 
at the resuit l|3.18|l . we have used two-tori which should lead to a conservative bound on t l . 
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of unit is given by v in Eq. ()3.18j) . both on the original space and on the mirror, the parameters 
are quantized in precisely the SciiïiG Wciy as c a cinci fi a j thcit is, 




(3.22) 



where ëj are integers. 

Finally let us comment on the other parameters which will appear in our discussion and that we 
did not discuss here. Given the above conventions ail the flux/torsion parameters are quantized in 
terms of the same unit, and we shall assume the same for the flux 6 and torsion parameters of the 
more gênerai models considered in section [2~21 This is far from being a rigorous treatment, but the 
most natural and straightforward assumption one can make in the absence of detailed knowledge 
about thèse manifolds. 



4 Vacua of the basic models 

In this section we study moduli stabilization for the four-dimensional model based on half-flat 
mirror manifolds, as introduced in section l2~Tl For clarity we start with a simplified version where 
we consider only one size modulus, T, and one shape modulus, Z, together with the axio-dilaton, S. 
Later on in this section we will generalize our discussion to arbitrary numbers of T and Z moduli. 
Throughout, we will focus on supersymmetric solutions of the above Systems. 

4.1 The STZ model 

For the simple three-field model with one Kàhler modulus T = t + ir, one complex structure 
modulus Z = z + i( and the dilaton S = s + ier, the Kàhler potential (|3.4j) - (j3.6j) specializes to 

K = -ln(S + S) -31n(T + T) -3\n(Z + Z) , (4.1) 

where we have set dm = dm = 1. The flux/torsion superpotential ()3.1U|) now simply reads 

W = eT + eZ + yZ 2 , (4.2) 

and, including the gaugino condensate term, we have 

W = W + ke~ cS . (4.3) 

The F-equations for this model become 

F T = e -^W = 0, (4.4) 

F z = e + ifiZ-^-W = 0, (4.5) 

F s = -kœ~ cS - —W = . (4.6) 

The solution to (|4.4f) implies that 

W = —t , (4.7) 

6 The quantization of flux parameters in the generalised half-flat models can be discussed in more détail by studying 
their third cohomology and homology. It is likely to be more subtle than assumed in this paper. 
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which is a real quantity. Inserting this into (|4.5|) we find 



ez = et , 

(4-8) 

fj,z = U . 

Recall that our model is valid only in the régime of large volume and complex structure and, in 
particular, we have t, z ^ 0. Therefore, the second of équations ()4.8j) implies the vanishing of the 
magnetic flux term, that is \i = 0. Let us absorb the constant c in the gaugino condensate potential 
by defining the quantifies 

x = es , y = ca . (4-9) 
Then, using (|4.7|) . Eq. Q4.6JI can be written as 

-2fcre-oo8y = ^ l , , 

3 (4.10) 

2kxe~ x svny = . 

The value x = is unacceptable, as it would correspond to the strong (gauge) coupling limit. 
Consequently we have to impose siny = which fixes y to y = nir for some integer n. Finally, 
calculating directly the real and imaginary parts of W in Eq. Q4.3|) by inserting (|4.8|) . y = nir 
and u = 0, we can evaluate the constraint (|4.7|) . Combining ail results we find the most gênerai 
supersymmetric solution of our model to be 



er = — eÇ 



1 3k „-i/4 



et = ez = (-l) n+1 ^e 
4 

x = 1, (4.11) 
y = nir . 

Let us discuss this resuit. Thèse équations fix t, z, s and a, and the same holds for er + e£, while 
the orthogonal combination remains a fiât direction. It is clear from the above expressions that, in 
order to be in the large radius and complex structure limits, the torsion and flux parameters e and 
e should be sufnciently small. However, as those parameters are quantized, the best we can do is to 
stick to their minimal, non- vanishing, values which corresponds to \ë\ = \ë\ = 1 in équations (|3.2Uj) . 
Even for this choice, we need a value of k bigger than 5 to arrive at t > 1 and z > 1. In other 
words, it is difncult to stabilize fields in the large radius and large complex structure région and, 
only by going to the limit of what one would consider reasonable parameter choices, can marginally 
consistent solutions be obtained. 

There is a similar problem with the gauge coupling since x is fixed at a relatively small value 
for the above solutions. Even using the relatively large E% beta-functions coefficient (|3.14j) we find 
for the inverse gauge coupling 

s = x /c= 1.19 , (4.12) 

47T 

which is barely in the weak coupling limit. 

7 Note that had we taken no gaugino condensate, that is k = 0, the above System admits a solution only if 
e = e = n = 0. This is the limit of compactifying on a normal Calabi-Yau manifold with no fiuxes turned on and it 
is in agreement with the resuit derived in Réf. (42) that the internai manifold has to be complex in order to obtain 
supersymmetric solutions. 
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4.2 The gênerai case 

Let us now briefly discuss the gênerai case, where h^ 1 ' 1 ' and M 2,1 ) are arbitrary integers. With the 
Kàhler potential as in Eqs (|3.4[) - (|3.6[) and the superpotential (|3.12[) . (|3.1U|) . we dérive the following 
F-equations 

Fjh = ei + KiW = 0, (4.13) 
F Z a = e a + i f i b d abc Z c + K a W = 0, (4.14) 

F s = -kœ~ cS - —W = , (4.15) 
2s 



with Ki and K a given by 



Ki = ~^d ijk tH k , K a = -Xd abc z h z c . 



Note that Ki and K a and, hence, W are real with the latter given by 

W = ^f. (4.17) 
As a conséquence, taking the imaginary part of Eq. (|4.14|) gives 

d ab cZ b V c = . (4.18) 

The matrix JC a b = d a b c z c is non-singular for a physical point z a in moduli space (as otherwise the 
Kàhler metric K ab would be singular at this point), so ail magnetic fluxes must vanish in order to 
have a supersymmetric solution. Equation l)4.15j) reproduces similar results to the case with only 
one T and one Z, namely y is constrained to take the values y = mr with n integer, while x must 
obey 

(-l) n+1 2kxe- x = ^f- . (4.19) 

As before, we can compute the value of the superpotential directly by inserting \x a = 0, y = nir 
and Eq. (|4.14j) into Eq. (j3.1flj) . On the other hand, we know from Eq. (|4.17j) that the imaginary 
part of W must vanish which leads to the constraint 

e l r i + e a C = . (4.20) 

This will be the only relation involving the axions, so we can only fix one of them while we are left 
with M 1 ' 1 ) + - 1 fiât axion directions. Matching the real part of W with Eq. (|4.17f) fixes the 
value of the dilaton to 

x = \ , (4.21) 

while the t % and z a moduli obey 

e . = ( _ ) n+l^ e -l/4^i|^ ^ 6a = (_)»+l^ e -l/4 ^ / , (4.22) 

It appears that generic analytic solutions to thèse last équations for t l and z a cannot be written 
down but, of course, solutions can be obtained, either analytically in simple cases or numerically, 
once explicit sets of intersection numbers dijk and d a b c have been fixed in the context of a particular 
model. We will not carry this out explicitly, as we have already seen that there exist fiât axion 
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directions and that the value of the dilaton is unchanged from the simple three-field case, so that 
weak gauge coupling is difficult to achieve. However, it is clear that solutions to Eqs (|4.22j) will be 
of the form t l ~ /c/(flux or torsion) (and similarly for z a ) so that flux/torsion quantization makes 
it hard to obtain vacua in the large radius and large complex structure limits. In summary, we 
have seen that the gênerai model shows ail the major problems that we have already found in the 
simple three-field case. 

Let us consider if there are any alternative ways around the above problems. Clearly, some of 
the difficulties arise because the supersymmetry condition forces us to set the magnetic fluxes ji a 
to zéro. This problem may not arise for non-super symmetric vacua. However, we note that the 
scalar potential only dépends on the combination e^r 1 of the axions r* (since this is true for the 
superpotential and the Kàhler potential is axion independent). Hence, even for non-supersymmetric 
vacua we will have at least h > — 1 fiât directions. A possible way forward could then be to study 
non-supersymmetric solutions for models with only one T modulus, that is, h 1 ' 1 = 1. We will not 
do this in the présent paper, as we focus on supersymmetric vacua, but we simply note that, for 
models with h ' = 1, there is still a chance for consistent (non-supersymmetric) vacua with ail 
moduli fixed. 

Another possibility is to modify the superpotential (|3.10j) to 

W = e{T l + idyfcm^'T* + e a Z a + % -d ahci x a Z h Z c , (4.23) 

that is, by including magnetic torsion terms with torsion parameters m 1 . Although this is a sug- 
gestive extension of the basic model, with a superpotential perfectly symmetric between the Kàhler 
and complex structure moduli parts, we do not currently know of a convincing dérivation of such 
a model in the context of the heterotic string. Given this situation, we will only give a very brief 
summary of the results for moduli stabilization we have obtained for such models. We find that 
there exist supersymmetric vacua with ail moduli stabilised and values of the dilaton x in the range 
x € [0, 1]. For suitable choices of parameters x ~ 1 can be achieved and, with the E$ beta-function 
coefficient (|3.14|) . this implies an inverse gauge coupling of at most s ~ 4.8. This is in the weak 
coupling région, although still well away from the "phenomenological" value s ~ 24. The values of 
t % and z a are proportional to the magnetic torsion/flux parameters, that is, t l ~ m 1 and z a ~ fi a , 
but the constant of proportionality in thèse relations is such that large radius and large complex 
structure can barely be achieved by minimal flux/torsion parameters and a value of k at the upper 
end of the reasonable range. In summary, adding a magnetic torsion term can solve two of the 
three problems of the basic model, namely fix ail moduli and generate weak coupling (although 
perhaps not to the desired extent), but achieving large radius and large complex structure remains 
problematic. 

Why is it so difficult to generate sufficiently large values of t 1 and z a ? In ail examples the 
values of thèse fields were basically determined by an expression of the form k / (flux or torsion) . 
The lower bound on the flux/torsion parameters due to quantization, combined with the fact that 
k is expected to be not too large in a' units, rules out large field values. The proportionality 
of the field values to the constant k in the gaugino condensate potential can be traced back to 
the fact that the flux/torsion part W of the superpotential does not have non-trivial globally 
supersymmetric solutions by itself. For example, from the superpotential (|3.10f) for the basic model 
we have W^i = which has no solution unless the torsion parameters ei vanish. On the other 
hand, if W had globally supersymmetric solutions, the values of t 1 and z a at this global level would 
be determined by fluxes and torsion parameters only. Provided the locally supersymmetric solution 
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can be obtained as a perturbation of the global one this would essentially découplé the values of f 
and z a from k and potentially solve our problem. To understand the local solution as a perturbation 
of the global one, the value Wo of the superpotential taken at the global solution should be small. 
A small | Wo| also facilitâtes weak coupling as should be intuitively clear from the structure of the 
superpotential (|3.12j) . We will explain those statements in more détail in the following section, 
where we analyze the models based on generalised half-flat spaces. As we will see, for those models 
W has, in gênerai, global supersymmetric solutions and, under certain conditions, |Wo| can indeed 
be made small. 

5 Generalised half-flat models 

In this section, we will analyze moduli stabilization for the generalised half-flat models introduced in 
Section 12.21 They are significantly more complicated than the basic models of the previous section 
as they involve more flux/torsion parameters per field. It will therefore be harder to find simple 
analytic solutions and we will have to use approximations and numerical methods. Also, the main 
part of our discussion will be within the simplest three-field STZ model, where M 1 ' 1 ) = M 2 ' 1 ) = 1. 
This model already contains eight flux and torsion parameters. However, our main results should 
carry over to the gênerai case, which we will discuss at the end of the section. 

5.1 Relation between locally and globally supersymmetric solutions 

Before we launch into the analysis of the STZ model, we would like to understand the relation 
between globally and locally supersymmetric solutions of our models in gênerai. This will also 
provide us with a practical way of finding supersymmetric vacua. We start with a superpotential 
of the form 

W = W + ke~ cS , (5.1) 

where W = W(T* , Z a ) is independent of the dilaton and stands for the flux/torsion part of W. 
For the purpose of the présent discussion we can keep W arbitrary but, of course, we have in mind 
the concrète form (|3.7|) . Let us assume that T l = Tq and Z a = Zq is a globally supersymmetric 
minimum, that is, it satisfies d{W(To, Zq) = d a W(To, Zq) = 0, and let Wo be the value of the 
superpotential at this minimum 

W = W(T ,Z ). (5.2) 

Further, let Mq be the typical moduli mass at this minimum, computed at the global level from 
the second derivative of W and let us assume that |Wo/-Mo| « 1. It is not hard to see that this 
condition is sufficient to ensure that the F-equations, Fi = and F a = 0, are approximately satisfied 
by the global solutions Tq and Zq, up to small corrections 

81* ~sz a ~^ . (5.3) 
M v ; 

Note that, in the above analysis, we have used the fact that we are working in the large radius 
and large complex structure limits, that is, the moduli fields t % and z a are bigger than one. Values 
much larger than one for thèse fields will make the approximation even better. For smaller field 
values (for example near conifold points) the above argument would have to be refined. 

The spécifie flux/torsion superpotentials and their derivatives will generically be of order one 
or larger due to the quantization of flux and torsion parameters. To satisfy the above condition it 
will, therefore, be sufficient to have |Wo| -C 1 in a' units. 
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Let us turn now to the dilaton F-equation 



F s = d s W + K S W = . (5.4) 

Expanding W around the global solution and using (|5.3|) . the leading contribution to the above 
équation will corne from 

F s ~ -cke~ cS - ^- (W + ke~ cS ) = , (5.5) 



which then yields the solution 



(2x + 1)6 



y = — arg 



2s 

W 



(5.6) 



for the rescaled dilaton 



x = es , y = ca . (5-7) 



It is this class of super symmetric solutions which we will be looking for in our generalised models. 
For such vacua the values of T l and Z a are determined at the global level, thereby potentially 
circumventing the problems of achieving large radius and large complex structure encountered in 
the previous section. Furthermore, it is clear from Eq. (|5.6|) that a small value |Wo| facilitâtes a 
large value of the dilaton x and, hence, weak coupling. 

Note that this procédure is slightly différent to the one outlined in Refs |HJ HUj . where the 
issue of integrating out heavy fields in SUSY théories was addressed and applied to the KKLT [2] 
scénario. Whereas in most papers the F-equations are used to integrate out heavy moduli we, as 
indicated above, start with a globally super symmetric solution and, by making |Wo| small, ensure 
that an approximate solution to the full F-equations exists. In fact, we have numerically checked 
our procédure and verified that a solutions of the full F-equations indeed exist close to the globally 
supersymmetric ones, provided |Wo| is small. An explicit example will be presented in the next 
subsection. 



5.2 The STZ model 

In this subsection, we discuss the three-field model with a single Kàhler modulus T, a single complex 
structure modulus Z and the dilaton S. From Eqs ()3.4[) - ()3.6|) the Kàhler potential reads 

K = -ln (5 + 5) -31n(r + f) -31n(Z + Z) . (5.8) 

The gênerai flux/torsion superpotential 1)3.7(1 becomes 

W = i(£ + ieT) + (e + ipT)Z + + iqT)Z 2 + \{p + irT)Z 3 , (5.9) 

2 o 

where we have chosen the signs of the flux parameters £, e, fi and p and the torsion parameters e, 
p, q and r for convenience. We recall those parameters are subject to a number of constraints (|3.8[) 
and (|3.9j) . However, the first set of thèse constraints (|3.8j) is trivially satisfied for h^ 1 '^' = 1 while 
the second set reduces to the single condition 

Çr — eq + pp — pe = . (5.10) 
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We remind the reader that this constraint originates from the relation dH = 0, which is a con- 
séquence of the Bianchi identity (|2.25j) if the order a! terms on the RHS cancel. This happens, 
for example, if the standard embedding is chosen, that is, if the gauge connection is set equal to 
the spin connection. On the other hand, for more gênerai gauge bundles the constraint can be 
avoided. We will discuss both cases, with and without the constraint (|5.1U|) . As usual, we take the 
full superpotential to be 

W = W + ke~ cS , (5.11) 

with W as in Eq. (j5.9|) . 

Following the procédure outlined at the beginning of the section, we will start by searching for 
global supersymmetric vacua of W. Thèse can be found from 

W T = - e + ipZ - \z 2 + i r -Z z = , 

1 5 1 (5.12) 

W z =e + ipT + (ip, - qT)Z + - (p + irT)Z 2 = . 

For r ^ 0, the first of thèse équations is a cubic in Z which can be explicitly solved using Cardano's 
formula. One finds that for each choice of the flux/torsion parameters there is exactly one solution 
with z > if and only if the discriminant of the cubic is positive. The second équation can then 
be solved for T in terms of Z. 

For r = the solutions to the previous équations take the simple form 
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2 



For a given set of parameters we can now compute Wo and check whether it is much smaller than 
one. However, the size of the flux parameter space is such that it is virtually impossible to carry 
out an analytic search of favoured régions. Therefore we resort to a numerical scan, varying the 
flux/torsion values (which are integers) in a certain range from —M, . . . , M. 

We have found the corresponding vacuum solution for each set of parameters, keeping only 
those vacua in the large radius and large complex structure limits, that is, with t > 1 and z > 1. 
We first carried out this procédure imposing the constraint (|5.10[) . For the case r^Owe find, for 
flux/torsion parameters in the range —20, . . . ,20, that |Wo| > 0.8 always. A similar lower bound 
in |Wo| is found for r = where we have searched in the range —70, . . . , 70. Furthermore, the 
lower bound for |Wo| is reached for relatively small values for the flux/torsion parameters and |Wo| 
does not decrease any further as the parameters range is increased. We take thèse results as strong 
évidence that small values of |Wo| cannot be obtained within this model and, hence, that it will be 
difficult to achieve weak coupling. 

Next, we have repeated the above procédure but without imposing the constraint (|5.1()j) focusing, 
for simplicity, on the case r = 0. The results are surprisingly différent from the ones obtained when 
the constraint is imposed. In particular, we find that vacua with small values of |Wo| can be found 
without any problem once the constraint is dropped. A useful way of summarizing the resuit is 
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to introduce the quantity N = N(M,w), defined as the number of vacua (with t > 1 and z > 1) 
found in the range —M, . . . ,M of flux/torsion parameters, with associated superpotential values 
Wo satisfying |Wo| < w. Numerically, we find that N is well-described by the scaling law 

N(M,w) = N M x w v , (5.14) 

where Nq, x and y are real constants. From a numerical search of ail flux/torsion parameters with 
M = 10, 20, ... , 70 we find that 

Aq-0.17, x~5.0, y -2.1. (5.15) 

This can be easily seen from the following two figures. In Figure^we have plotted the total number 




log(M) 

Figure 1: Total number of vacua, N, as a function of the range of flux/torsion parameters, M (in 
logarithmic units). 

of solutions with \Wq\ < 1 (i.e. taking w = 1 in équation ([5.14)0 for M = 10, 20, . . . , 70. In Figure[2] 
we have presented the resuit for N(70,w), the number of vacua in the flux/torsion range from 
—70, . . . , 70, as a function of the superpotential value, \W\. The numerical value of x can be easily 
understood intuitively. Although we consider a model with seven parameters, the requirement of 
small Re(Wo) and small Im(Wo) effectively fixes two of thèse parameters, leading to a scaling law 
with power 5. We also remark that the value of y close to y = 2 corresponds to a nearly uniform 
random distribution of Wo values in the complex Wo plane. 8 

As can be seen from Figure |Wo| values of 0.01 and smaller can be obtained. Our resuit 
gives an indication of what fraction of vacua leads to a gauge coupling of s ~ 24 as suggested by 
gauge unification in the MSSM. Assuming the E$ beta-function coefficient (|3.14j) such a value for 
s translates into x ~ 5 which, from Eq. ()5.6|) (setting k = 1 for simplicity) implies |Wo| — 1/14 
(or, equivalently, log|Wo| — —1.15). This means approximately a fraction of 10 -3 of ail vacua with 
|Wo| < 1 lead a gauge coupling sufficiently weak to be compatible with gauge unification. For a 

8 We thank Nuno Antunes for pointing this out to us. 
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2 1 - 

-3 



log(|W|) 



Figure 2: Number of vacua N with | Wo| < w as a function of w (in logarithmic units), for 
flux/torsion parameters in the range —70, . . . , 70. 



condensing gauge group smaller than Eg or a value of k smaller than one this fraction will decrease 
accordingly. 

We have also checked the case r ^ without the constraint Q5.1UJI and the results are similar 
to the one obtained for r = 0. 

Finally we would like to show a numerical proof of the validity of our procédure (i.e. the use 
of the global SUS Y condition W;(T , Z ) = W a (T , Z ) = with |W /M | < 1 in order to find an 
approximate solution for the T and Z fields). We have chosen, within this r = case, values for 
the flux/torsion parameters as follows 



e = -7, e = -4, fi = 2 
for which Eqs (j5.13|) give the field values 

Zo 
to 
Co 

70 



P = 5 , p=l , q = 2 , £ 



2.598 , 
2.165 , 
0.500 , 
2.250 , 



13 



(5.16) 



(5.17) 



with |Wo| = 0.167. As it is illustrated in Figure 3, thèse values are very close to the actual solution 
to the F-equations, given by 



to 
Co 

70 



2.598 , 
2.164 , 
0.520 , 
2.300 . 



(5.18) 
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Using the beta-function coefficient Q3.14|) and k = 2, the values for the dilaton field are found to 
be 

s =23.190 , 

5.19 

a =22.399 , 

which proves that phenomenologically viable values for the gauge coupling can be obtained for 
reasonable values of the flux parameters. We should add that we have computed the Hessian 





Figure 3: Contour plot of the potential, on the (t, z) plane, for the example shown in the text, 
Eq. (jQH). The solid (dashed) lines are the Re(F T ) = Re(F z ) = (Re(Wr) = Re(W z ) = 0) local 
(global) SUSY conditions. Note that the lines Re(Wr) = and Re(-Fr) = coincide. 



matrix for this example and we have explicitly checked that ail eigenvalues are positive at the point 
where the F-equations vanish. This is quite important as the potential at the minimum is négative, 
actually given by V = — 3e^|VVo| 2 , and very small, of the order of 10 -7 . This means that, in a 
small région around the true minimum, the potential will shift from négative to positive values 
and, this being a mufti-variable potential, it is easy to get mistaken about the real position of the 
minimum. For example, if we were to use the global supersymmetric solution given by (j5.17jh due 
to its closeness to the real one, (|5.18|) . the no-scale cancellation mechanism would take place and 
the scalar potential would read, at this point, V = +3e^|Wo| 2 . That is, we would predict a dS 
vacuum (of order 10 -7 ) were, in reality, the only minimum in that région is AdS. 

Before we conclude this section we would like to make a few comments about the gênerai case, 
with an arbitrary numbers of moduli fields. As it is évident from (|2.3U|) and ()2.27j) . the number of 
parameters grows rapidly with the number of fields, making a numerical search for vacua infeasible. 
We have, therefore, not attempted to extend the numerical analysis beyond the three-field STZ 
model. However, the expérience with this model shows that dropping the constraint (|3.9j) is crucial 
in order to obtain small values of | Wq| and we expect this to be true in gênerai. Conversely, 
multi-field models without the constraint (|3.9j) should be at least as flexible as their three-field 
counterpart and should, therefore, allow small |Wo| values without problems. From our gênerai 
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argument relating globally and locally supersymmetric vacua, this should then allow consistent, 
locally supersymmetric vacua at weak coupling as determined by Eq. Q5.6|) . As for the scaling 
law (|5.14|) . in the gênerai case one would expect scaling powers x ~ n — 2, where n is the total 
number of parameters in the model and y ~ 2, leading to the same uniform random distribution 
in the W plane which we have observed in the STZ model. 

6 Conclusion 

In this paper we have analyzed the vacuum properties of various classes of heterotic models on 
certain manifolds with SU(3) structure. After a review of the heterotic string on half-flat mirror 
manifolds [32], defined by (|2.2U|) . we have derived the superpotential for a more gênerai class of 
manifolds with SU(3) structure which were introduced in Refs |34| I35| VAty . We have explicitly 
verified in thèse models that the application of the heterotic Gukov-Vafa-Witten type formula for 
the superpotential leads to the same resuit as an explicit réduction of the ten-dimensional bosonic 
terms. The resulting superpotential, which is given in Eqs (|2.30|) . ()2.27j) . resembles very much 
the one obtained in type IIB orientifold compactifications suggesting that one may recover the 
nexibhity of type II models in the heterotic case. Thèse flux/torsion superpotentials dépend on 
h^ 1 ' 1 ' Kàhler moduli T % as well as on M 2 ' 1 ) complex structure moduli Z a , but are independent of the 
dilaton S. We have, therefore, supplemented our superpotential with a contribution from hidden 
sector gaugino condensation in order to stabilize the dilaton, which is again similar to the type IIB 
constructions where non-perturbative terms need to be added in order to fix the Kàhler moduli. 

We have first analyzed moduli stabilization for the models based on half-flat mirror manifolds 
and have found a number of problems. Generally in those models, M 1 ' 1 ) + M 2 ' 1 ) — 1 axion directions 
remain flat, and it is hard to achieve the large radius and large complex structure limits as well as 
weak gauge coupling. In models with additional magnetic torsion terms, the flat axion directions are 
lifted and moderately weak coupling can be achieved, while stabilizing field in the large radius and 
large complex structure limits remains a problem. However, such models with additional magnetic 
torsion terms, although a plausible extension of half-flat mirror models, cannot currently be derived 
within the heterotic string. 

We have traced the root cause of the aforementioned problems to the fact that the flux/torsion 
superpotential W is too simple to allow for globally supersymmetric vacua. Consequently, we have 
analyzed moduli stabilization for some generalised half-flat models whose associated superpotential 
is significantly more complicated. We have seen that consistent weak-coupling vacua can be ob- 
tained if the flux/torsion superpotential has globally supersymmetric vacua with a superpotential 
value Wo satisfying | Wo| <C 1. The value of the dilaton and, hence, the gauge coupling, is then di- 
rectly related to | Wo| ■ We have verified that the superpotential for the generalised half-flat models 
has indeed globally supersymmetric vacua with ail Kàhler and complex structure moduli stable. 
However, the requirement of small | Wo| turned out to be more subtle. For the standard embedding 
of the spin connection into the gauge group, the resulting Bianchi identity dH = for the NSNS 
form led to a constraint (j3.9j) on the flux/torsion parameters which ruled out the possibility of small 
Wo, at least within the range of flux/torsion parameters covered by our numerical scan. However, 
for more gênerai gauge bundles, the constraint should be dropped and vacua with small |Wo| can 
easily be obtained in this case. The number of such vacua as a function of |Wo| for a simple model 
with three fields, S, T and Z, has been plotted in Figure [2j Using Eq. (|5.6|) one can estimate 
that, typically, the fraction of vacua that lead to a sufficiently weak gauge coupling consistent with 
gauge unification, is 1CP 3 . Our results establish the existence of consistent, weak-coupling AdS 
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vacua within generalised heterotic half-flat models. 

In the liglit of thèse results, it is clearly désirable to get to a better understanding of half-flat 
compactifications and their gêner alizations, in particular with regard to the précise nature of the 
manifolds involved, the rules for quantizing flux and torsion parameters in those compactifications 
and the inclusion of gauge and gauge matter fields. We will leave those tasks for future publications. 

Acknowledgments We would like to thank Nuno Antunes, Kiwoon Choi, Jan Louis, Hans-Peter 
Nilles and Silvia Vaulà for helpful discussions. BdC, AL and AM are supported by PPARC. S G is 
supported by the JSPS under contract P03743. 

Appendix 

A Ricci scalar for the "extended half-flat" manifolds 

This section contains a generalization of the resuit obtained in the appendices of Réf. |HS]) where 
the Ricci scalar for half-fiat manifolds mirror to Calabi-Yau with NS-NS fluxes was computed. Here 
we will follow this calculation closely, by recalling the main identities which remain valid, while 
pointing out the places where it differs from the simple half-flat case. 

To set the stage, let us briefly recall a few features of manifolds with SU(3) structure. Such 
manifolds are characterized by the existence of an almost complex structure with the associated 
fundamental form J and a (3, 0) complex form Çî which are invariant under the action of the SU(3) 
structure group. More concretely this means that the forms J and fl are covariantly constant 
with respect to some connection V (T) , which in gênerai has a torsion. Decomposed into SU(3) 
représentations the torsion falls into five différent classes Wi , . . . , W5 which are given by 

dJ = Wi^ + W 4 AJ + W 3 , (A.l) 
dn = WiJA J + W 5 AÎÎ + W 2 . 

Since the torsion on manifolds with SU(3) structure measures the departure from Calabi-Yau 
manifolds (which are Ricci fiât) it is clear that the Ricci scalar of the SU(3) structure manifolds 
dépends on their torsion. Thus, in order to compute the Ricci scalar, we will need to know ail the 
components of the torsion and for this we will use équations (|A.1|) above and the relations 

dn = EiCj 1 , (A.2) 
dJ = (f PAi )[3 A - {tqt)a A , (A.3) 

which are easily derived from (|2.4j) and (|2.21|) . with the quantifies Ei defined in équation (|2.34|) . 
Note we have postulated that the basis forms in the above équations have the same SU(3) properties 
as their Calabi-Yau counterparts. Hence, a A and (3 A in Eq. (|2.21j) are primitive and, consequently, 
W4 has to vanish. Moreover as dVL in Eq. (|A.2j) is a (2, 2) form, W5 also vanishes. The other torsion 
components T\ + 2 and T3 are found to be 

(ïi+a)^ = 4p|2(û'Wn^7. ( A - 4 ) 

(T 3 ) a fô = -^t l (pAiPa^ ~ qfaAaPj) ■ (A.5) 
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Here ||f2|| 2 = gîî a/ g 7 f2 Q ^ 7 is a function of the complex structure moduli which is related to the 
Kàhler potential for thèse fîelds via V||f2|| 2 = e~ Ki ~ Z \ where V is the volume of the manifold. Note 
that, in this case, the quantities Ei are neither real nor constants as it happened for the case of 
half-flat manifolds and thus one generically has 

T G Wi © W 2 © W 3 • (A.6) 

However it is important to note that the nature of the indices of the torsion components is the same 
as in the half-flat case, and the torsion itself is still traceless. As a conséquence, the expression of 
the Ricci scalar in terms of the torsion obtained in |33| still holds 

R = (T 1+2 ) Q/37 (T 1+2 )^ - 6(T 1+2 )^ 7 (T 1+2 )^ Q + (T 3 ) a ^(T 3 ) a ^ + ce. (A.7) 

To obtain the integrated Ricci scalar we perform the same steps as in Réf. [33j and, using the 
relations (jB.lOj) . we obtain 

/ = ^(J'EiÈj-tflPEiÉj) 

Jx 

-\M(pAi - q?MoA)(lmM- 1 ) AB (PBj ~ <$Mbd) ■ (A.8) 

After taking into account various coefficients and rescalings, the contribution of gravity to the 
potential in Einstein's frame can be rewritten as 

e~ K V g = 8(g ij EiËj - 4tH j EiËj) (A.9) 
-Ae~ K(Z) [lm(ë A ) - lm(f)McA) (imAT 1 )^ [lm(ë B ) - Im(^ D )M DB ] ■ 



B Some useful results on spécial geometry 

As it is well known, the moduli space of Calabi-Yau manifolds splits into a product of two spécial 
Kàhler manifolds, one for the complexified Kàhler class déformations and one for the complex 
structure déformations. Since thèse geometries are at the heart of the four-dimensional physics 
obtained from compactifications on Calabi-Yau manifolds we review in this appendix some of the 
properties of the spécial Kàhler manifolds which we need in the main part of the paper. We mainly 
follow Réf. @B]. 

The main feature of spécial Kàhler manifolds is that their geometry is completely determined in 
terms of a holomorphic function 7i, called the pre-potential. In terms of the projective coordinates 
X p where P = 0, . . . ,n (n being the complex dimension of the manifold), the pre-potential is a 
homogeneous function of degree two, that is, it satisfies X p Tip = 27ï, where Hp = g^pH. In fact, 
one does not always need to rely on the pre-potential and it may be sufficient to work with the 
period vector 

° = (£)- < R1 > 

Let us further introduce the symplectic inner product <, > as 

<O,Ô>=O T (^ 1 ~^Ô = {H P X P -H P X P ) . (B.2) 

With this notation the Kàhler potential can be written as 

K = - log (i < O, Ô > ) = - log [i(H P X p - HpX p )] , (B.3) 
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while the Kàhler metric is given by the usual formula 



gp q = OpBgK . (B.4) 

Here the derivatives are with respect to the affine coordinates X p = X p /X°, where p,q,--- = 
1, . . . , n. It is also useful to introduce the Kàhler covariant derivative of the periods O 



e 



- K/2 U q = e~ K/2 ( ^j=V p O = {d p + d p K)0 . (B.5) 



The period matrix Q, which is a complex symmetric matrix is now required to satisfy 

(B.6) 

hpq =Qpsf q ■ 

It can be shown, see Ref [44] . that, in terms of the pre-potential 7i, the period matrix has the form 

Qpq - H PQ + 2z {lmH)RsXRx s • (B-7) 

where we have denoted Tïps = âf^âf 3 "^' 

With this one can prove the following relations 

<o,u p > = < O, Ûp >= , 

g p g = -i<U p ,Ûg>=-2f p p lmQ PQ f9 , (B.8) 

In order to make this discussion less abstract let us first apply the above formalism to the com- 
plex structure moduli space of Calabi-Yau manifolds. The periods O are now determined by the 
holomorphic (3, 0) form Çl. Moreover the inner product ()B.2|) becomes now the inner product for 
three-forms on the Calabi-Yau manifold. With this one immediately finds that the Kàhler potential 
(|B.3|) precisely reproduces the one from équation (|2.14|) . Finally, the Kàhler covariant derivatives in 



équation (|B.5|) give the components of the (2, 1) forms in the basis (|2.2j) . With thèse identifications 
it is easy to see that most of the relations in équation (|B.8|) are straightforward, the only non-trivial 
ones involving the period matrix. Denoting the period matrix by M in this case and the indices 
P,Q, ... by A, B, ... = 1, ... , h^ 2,l \ the last relations in équation (|B,8|) become 



g al = -2/flmAWl 



b ' 



1 



e K D a Z A D- b Z B = iïffg ab = -^ImM)" 1 AB - e K Z A Z B . (B.9) 
Finally we note that in the basis 1)2. 2JI the period matrix A4 can be found to be 
[ a A A*f3 B = -(ReM) A c(tmMy 1CB , 

[ a A A*a A = -{lmM) AB -(ReM) A c{lmM)- lCD (ReM) DB , (B.10) 
(3 A A *(3 B = -{JmM)~ 1AB ■ 



Y 3 



25 



While the pre-potential is typically a complicated function it simplifies considerably in certain 
limits in moduli space, such as large radius and large complex structure limits for the Kàhler and 
complex structure moduli spaces, respectively. In those limits, the pre-potential is given by a cubic 
function 

1 d var X p X p X r 

Such a cubic pre-potential defines what is known as very spécial geometry. Writing the affine 
coordinates as 

XP = e + ix p , (B.12) 

one finds for the Kàhler potential 

K = -]n(^]C\ , (B.13) 

where 

K = U pqr (X p - X p ){X q - X q ){X r - X r ) = d pqr x i x j x k . (B.14) 
8 

From équation ()B.7|) one can also define a period matrix in this case and then the relations 
()B.8|) follow by straightforward algebraic manipulations. There are a number of further very spécial 
geometry relations which are useful in the main part of this paper. First, let us define 

/C-£) — d'pq'p x x - ^-'pcj — drpqy x é ^ 13 .15^ 



With this, the first derivatives of the Kàhler potential K with respect to X p , denoted by K p and 
the Kàhler metric K pq can be written as 

K P = — jj^, (B.16) 

Defining fields x p = g pq x q with lowered indices it is easy to show that 

x p = ^, x pX v = \. (B.18) 

Thèse formulae lead immediately to the "no-scale" relation 

K pq K p Kg = 3 . (B.19) 



Finally we note that a typical flux superpotential, for example as it arises from the Gukov-Vafa- 
Witten formula, can be written as 

W = epX p — m p Hp . (B.20) 

Here ep and m p dépend on the fluxes and can be either real constants or can also dépend holo- 
morphically on other (super)fields in the theory, but not on X p , as we have seen in section l2~2l For 
the cubic pre-potential ()B.11|) the dependence on the physical degrees of freedom X p can be made 
explicit after setting X° to one and we find 

m® 

W = e + e p X p + -d pqr m p X q X r - —d pqr X p X q X r . (B.21) 

Note that after transforming to the "phenomenological" convention for X p by X p — > —iX p (and 
after dropping an overall factor of —i from W which is irrelevant) the constant and quadratic terms 
in the above superpotential pick up a factor of i. 
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C Supergravity conventions in d = 4 and stability of supersym- 
metric vacua 



In this appendix we summarize conventions and relevant formulae for four-dimensional N = 1 
supergravity |45j . Further, we présent an elementary proof that solutions to the F-equations are 
always stable vacua. 

The bosonic terms in the action of four-dimensional N = 1 supergravity coupled to chiral fields 
= (S,T\Z a ) read 



S = 5- / \J — g d x 



R + K XY d^ x d^ Y + v 



(Cl) 



where K4 is the four-dimensional Newton constant. As usual, Kxy = dxdyK is the Kàhler metric 
while the potential V is given by the standard formula 

V = e K (K XY F X F Y - 3\W\ 2 ) , (C.2) 

where K XY is the inverse Kàhler metric and the F-terms Fx are defined by 

F X = W X + K X W . (C.3) 

Here a subscript X dénotes a derivative with respect to 4> x , as usual. Note that we have considered 
the chiral fields and the Kàhler potential to be dimensionless, while the superpotential has dimension 
one, a convention which is convenient for the discussion of moduli fields and in line with the formulae 
in the main part of the paper. 



In this paper we were interested in supersymmetric vacua of the potential (|C2|) . that is, vacua 
which can be obtained by solving the F-equations 

F x = W X + K X W = . (C.4) 

It is easy to show, from Eq. (jC2|) . that solutions to the F-equations indeed constitute extremal 
points of the potential V. The cosmological constant, Vq, at such an extremal point is given by 

V = -3e K \W\ 2 . (C.5) 

Without fine-tuning (to make W at the extremal point vanish) this value will usually be négative 
and, hence, we are generically dealing with AdS vacua. The stability of AdS vacua in gravity 
coupled to scalar fields was analyzed a long time ago |461 147j by Breitenlohner and Freedman and, 
independently, by Abbott and Deser. They found that such vacua are stable if ail scalar field 
masses are larger than a certain lower bound, which is basically given by the cosmological constant 
Vq. Hence in AdS space, unlike in Minkowski space, négative (square) masses do not necessarily 
indicate an instability. In fact, it can be shown in gênerai, that this bound is always satisfied 
for supersymmetric vacua of supergravity théories and, hence, such vacua are always stable. We 
will now présent an elementary proof of this statement. 

Let us first formulate the Breitenlohner- Freedman bound for a theory with canonically normal- 
ized real scalars x l an d a potential V = V(x 1 )- We assume the potential has a stationary point at 
X* = Xo w hh négative cosmological constant Vo = V(xh) < 0- Define the mass matrix as 
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According to Breitenlohner and Freedman this stationary point leads to a stable, AdS vacuum if 
o-i > 3Vo/2, where Oj are the eigenvalues of M. A sufncient criterion for this to be the case is that 

£ T MÇ - ^ Ç T e > , (C.7) 

for ail vectors £ in field space. To see that this inequality implies the one for the eigenvalues, choose 
£ to be the eigenvectors of the mass matrix M. It is useful for the application to supergravity to 
re- write this criterion for non-canonical kinetic terms Gijd^d^X 3 1 where Gij is the metric on field 
space. Then, Eq. (jC.7|) takes the form 

Ç T MÇ - ^ Ç r G£ > , (C.8) 
where the mass matrix M is now, of course, defined with respect to the non-canonical fields. 



We would like to apply the criterion (|C.8|) to the case of four-dimensional N = 1 supergravity 
with complex scalars (j) X , Kàhler potential K and superpotential W . We consider a solution <p x of 
the F-equations (|C4|) . This solution is automatically a stationary point of the potential V and it 
préserves supersymmetry. The cosmological constant Vq at such a vacuum is given by Eq. ()C.5|) . 
From Eq. (RT2|) one finds for the second derivatives of V at Fx = after a bit of computation 



Vxy\ f=0 = -e K WF XY , (C.9) 
K zT F zx F f y-2K X y\W\ 2 ] , (CIO) 



y - I P K 
v XY\F=0 ~ e 



where Fxy is the derivative of Fx with respect to <fr . Combining thèse results (and taking care 
to convert real into complex expressions) it is then straightforward to compute the LHS of the 
criterion ()(J.8|) which takes the form 



2K zf F zx F TY i x ê - WF XY i X f ~ WF XY i X ê + \\W\ 2 K X H X ^ 



= e K (-L\VK ZX Ç X - y/2F zx çA K zf LLwK fY f - y/2F?yçA. (Cil) 

The last line is obviously positive and, hence, the criterion is satisfied. The conclusion is that 
any supersymmetric AdS vacuum of four-dimensional N = 1 supersymmetry satisfies the Breiten- 
lohner /Freedman criterion and, therefore, constitutes a stable AdS vacuum. 
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